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We propose an optical ring interferometer to observe environment-induced spatial decoherence of 
massive objects. The object is held in a harmonic trap and scatters light between degenerate modes 
of a ring cavity. The output signal of the interferometer permits to monitor the spatial width of the 
object's wave function. It shows oscillations that arise from coherences between energy eigenstates 
and that reveal the difference between pure spatial decoherence and that coinciding with energy 
transfer and heating. Our method is designed to work with a wide variety of masses, ranging from 
the atomic scale to nano-fabricated structures. We give a thorough discussion of its experimental 
feasibility. 



I. INTRODUCTION 

The Schrodinger cat is a well known example of the dif- 
ficulty in clearly defining the border between the classical 
and quantum worlds. In this example, quantum theory 
allows for a macroscopic superposition of a dead and a 
live cat to exist while we have never been able to observe 
such a macroscopic superposition in nature. Indeed, this 
enigma has been the source of a century long debate. 

To the best of our knowledge, only few attempts have 
been made so far to artificially create macroscopic su- 
perposition states. One attempt dealt with a superposi- 
tion of two states of a multi-photon cavity field p|. A 
second dealt with a superposition of a magnetic flux di- 
rection, formed by two macroscopic counter propagating 
electron currents Q- Spatially separated superpositions 
of a trapped ion have been prepared and probed Q , and 
Bose Einstein condensates have been examined Q . Work 
on handedness of chiral molecules may also be considered 
relevant to this topic (5J- Microscopic, mechanical oscil- 
lators have been extensively discussed as well [j| QJ , and 
experiments are now entering the regime where quantum 
effects become observable y|. 

In this paper, we discuss a general oscillating massive 
object and its decoherence in the position basis, the basis 
which stands at the base of our classical perception. This 
spatial decoherence, also called "localization", is pre- 
dicted by numerous models for environment-induced de- 
coherence that are put forward to explain the appearance 
of classical reality from an underlying quantum world [9j : 
objects become localized in position due to their interac- 
tion with the environment. Of special interest is "pure" 
decoherence where localization can happen even without 
the transfer of energy, e.g., in a double-well potential. 
In order to reproduce the absence of macroscopic super- 
positions, most decoherence models use the mass of the 
decohering object as a central parameter along with pa- 
rameters such as time and spatial splitting of the super- 
position. Indeed, recent seminal experiments on matter 



wave diffraction have tried for the first time to explore a 
region of mass values beyond the usually experimentally 
feasible masses of elementary particles [l0| . The visibil- 
ity loss in the diffraction pattern is then a measurable 
signal of spatial decoherence. 

Interference experiments with freely propagating ob- 
jects, however, become increasingly hard to perform with 
larger masses for two main reasons: first, the de Broglie 
wave length becomes smaller and the required diffrac- 
tion gratings become more difficult to fabricate. Second, 
the spatial superposition created by diffraction is increas- 
ingly sensitive to decoherence because it is more difficult 
to isolate the freely propagating system from its environ- 
ment. This renders a controlled experiment with large 
masses extremely difficult. 

In this paper we propose an interference experiment in 
which no grating and no free propagation are needed. In 
fact, we avoid all together the need to create a well sep- 
arated spatial superposition, and hence the way should 
be open to perform controlled decoherence experiments 
with large masses. Our experiment is based on an optical 
interferometer to probe the state of an oscillating mass, 
e.g. a nano bead or a mechanical oscillator. Cont rary 
to previous work concerning oscillating mirrors 0, 1 1 1| - 
a symmetric ring interferometer is used, and the oscilla- 
tor is not required to have a high reflectivity. This fea- 
ture is advantageous when very light (thin) nano objects 
are investigated, as high transmittance does not pose a 
problem. We show that this setup can distinguish be- 
tween different models of decoherence dynamics so that 
information about this subtle process can be obtained 
experimentally. 

The next section describes the experimental setup. 
The theoretical model is presented in Sect. Ill, and solved 
approximately in Sect. IV. In Sect. V the experimental 
requirements are evaluated for specific examples. Gen- 
eral conclusions are put in Sect. VI. 
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II. EXPERIMENTAL SETUP 
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FIG. 1: Setup of the experiment. An atom or a more massive 
oscillating object, e.g. a nanoparticle, is held in a harmonic 
potential (P). An incoming photon is split and "hits" the ob- 
ject from both sides. The source (S), the beam splitter (BS) 
and the phase shifter (PS) form a preparation system which 
ensures that the photon mode is symmetric with respect to 
the symmetry axis. (PS compensates the phase difference be- 
tween the reflected and transmitted wave at the beam split- 
ter.) The same system acts as a detection system whereby 
the antisymmetric photon mode is sent to detector D. The 
potential is symmetric about the symmetry axis as well. If, 
for example, the object is initially prepared in a state of well- 
defined parity (e.g. the ground state of the potential) , its final 
state will remain a parity eigenstate unless decoherence breaks 
the initial symmetry of the photon+object system. 



The experimental setup we propose is sketched in 
Fig. ^ an object of mass m is confined in a potential 
(P) which can for all practical purposes be arranged in 
such a way that only the motion along one direction is 
relevant, the x-axis, say. We assume in this paper a har- 
monic confinement, V(x) — mil 2 x 2 /2. The object is held 
at the center of a ring interferometer, formed by massive 
mirrors (M) and an in-out coupling mirror (CM). The 
symmetry axis of the harmonic potential coincides with 
the plane of the beam splitter (BS), and hence with the 
symmetry axis of the whole setup. The BS and phase 
shifter (PS) prepare a symmetric photon mode, which 
excites a superposition (with zero phase shift) of right 
and left circulating travelling waves in the ring cavity. 
The BS and PS also act as a measurement apparatus for 
the outgoing modes, sending all anti-symmetric modes 
into the detector D. This is so because anti-symmetric 
modes of the cavity, are constructed from right and left 
circulating beams with a 7r phase shift. These are cou- 



pled out of the cavity through the CM in two separate 
directions. The PS reduces their phase difference to 7r/2 
as in a normal Mach-Zchnder apparatus, whereby the 
BS determines constructive interference in the direction 
of D. 

The experimental procedure we have in mind is the 
following: at t < the object is prepared in some equi- 
librium state at a given temperature. At t = the prepa- 
ration stops, and the environment, be it a thermal bath 
or some tailored environment, becomes dominant in the 
dynamics of the system. At t = t p , a probe light pulse is 
sent into the ring interferometer and interacts with the 
object. The consequent measurement of the probe pulse 
which exits the interferometer after the interaction, de- 
termines if the even symmetry of the initial photon has 
been altered. In principle also the energy change of the 
outgoing photon may be measured, but in this paper we 
make no use of this option. 

We note that the operation of the ring interferometer is 
not qualitatively affected when working with a "transpar- 
ent" object like an atom or a weakly scattering nanopar- 
ticle. For the ease of demonstration we first describe a 
single atom. The potential could be provided in this case 
by a magnetic trap [l2j. Following that, we extend the 
treatment to a massive nanoparticle to show how the ex- 
perimental sensitivity scales with mass and temperature. 

To conclude, the scheme provides an experimental 
measure of decoherence which may be used for probing 
large mass objects. It is interesting to note that future 
technologies may enable confinement of massive objects 
also in potentials other than harmonic, in which case de- 
coherence can lead to a clearer signature of spatial lo- 
calization without energy transfer ( "pure" decoherence) . 
However, in this work, we show that already for the most 
feasible of large mass potentials (the harmonic well), the 
ring interferometer is able to distinguish between differ- 
ent decoherence scenarios, and hence we truly present 
a realizable scheme for the probing of decoherence with 
large masses. 



III. MODEL 

We use standard techniques of open system quantum 
dynamics to model the experiment sketched in the pre- 
vious section. The initial state of the system (atom or 
nanoparticle) is described by a density matrix pa, which 
represents thermal equilibrium in the harmonic trapping 
potential at temperature Tq. Up to some probing time 
t = t 

Neumann equation 



p , the system evolves according to a Liouville-von 



PA 



--[H s ,p A ] +C\pa], 



with the harmonic oscillator Hamiltonian 



Ha 



2m 



(1) 



(2) 
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and a dissipative functional C that describes the influ- 
ence of the environment on the system. Its expression is 
detailed below. 

Around time t p , a light pulse is injected into the cav- 
ity to probe the state of the oscillator. The pulse sent 
into the cavity has a center frequency close to a cavity 
resonance and a narrow bandwidth compared to the free 
spectral range (FSR, which will be denoted by v). We 
assume that the oscillation frequency Q, is much smaller 
than v so that the interaction of the system with the 
field can only couple degenerate cavity modes. For sim- 
plicity, we also neglect the light scattered by the object 
into higher transverse cavity modes. The field in the cav- 
ity can then be described by only two degenerate modes 
with even and odd symmetry. As shown in Appendix A, 
the interaction Hamiltonian is given by 

Hint = hg ((aide - a' a ) cos(2fcr) + 

+ (a\a + a\a e ) sin(2fca;)) (3) 

where a e and a are the boson operators for the even 
and odd cavity modes, respectively, k = 2ir/\ is the cav- 
ity wave number, and g is the coupling strength. The 
excitation of the even mode is governed by the term 

H pump = -ih (alr)(t) - rf(t)a e ) , (4) 

where the pump amplitude rj is an operator whose state 
allows to describe any photon statistics of the incoming 
pulse. The leakage of photons out of the cavity is deter- 
mined by the transmissivity of the coupling mirror (CM). 
This process is accounted for by the Liouville functional 

C cav [p A F] = 2k ( a iPAFa\ - -[a\a u p AF ] + ) , (5) 

i— e,o ^ ' 

where k is the finite cavity linewidth (HWHM), [•,•]+ 
denotes the anti-commutator, and paf is the joint den- 
sity operator of the system and the cavity field modes. 
We assume for simplicity that all photons leaking out are 
actually detected with unit efficiency. Thus, the photon 
rate at the "odd" detector (D) is given by 2k times the 
number of odd photons in the cavity. 

As a signal, we will consider the average photon num- 
ber detected within in the time interval [t p ,tf], given by 



The functional 

C dec [pA]^~^[x,[x,p A }] (7) 

describes the limit of pure spatial decoherence (D is the 
classical momentum diffusion coefficient). It corresponds 
to the limit of weak friction and high-temperature envi- 
ronment in a Caldeira-Leggett model (see, e.g., [3). We 
obtain a pure thermalization Liouvillian by adding the re- 
quirement of detailed balance to the complete positivity. 
In terms of the creation and annihilation operators of our 
harmonic system Hamiltonian J2J), 

C therm [pA] = ^ _ ^b,p A ]^j (8) 

The stationary state of this functional is a canonical en- 
semble, with temperature ksTo = hi}/ ln(T i/T^) P.ll3|. 
In both cases, the description in terms of a master equa- 
tion is a reasonable choice if the coupling to the environ- 
ment is weak. 

Pure spatial decoherence can localize a system without 
the transfer of energy, for example in a deep double-well 
potential or through recoil-free scattering of probe par- 
ticles. In quantum Brownian motion, this corresponds 
to the limit where the environment correlation time is so 
short that the rotating wave approximation with respect 
to the system's oscillation frequency cannot be made and 
terms like b 2 and W 2 are retained in the interaction with 
the environment 0] . These terms drive the system from 
the initial ground state p A (0) — |0)(0| to a squeezed 
state. Therefore, after some time t p > 0, its density ma- 
trix in the energy basis will contain higher excited states 
such that off-diagonal elements become populated. Ther- 
malization (|SJ|, on the other hand, simply redistributes 
the weights of the diagonal elements (n\p A \n). The co- 
herences of the energy eigenstates produce a breathing 
motion of the spatial density which can be detected by 
the probe pulse as t p is varied. This signal then dis- 
tinguishes spatial decoherence from thermalization, and 
allows us to probe pure decoherence for a massive body. 



N = 2k Tr (ala oPAF (t)) dt . (6) 
Jt p 

The upper bound of the integration, tf , is chosen such 
that the detected photon rate is negligibly small by this 
time. The pulse length, and hence the detection window, 
is kept short compared to the decoherence timescale of 
the system. 

We focus in this paper on two different dissipative func- 
tionals that describe opposite extremes of dissipation and 
decoherence. We require the time-evolution to be a com- 
pletely positive semi-group, i.e. of Lindblad type |13| . 



IV. ANALYTICAL SOLUTION 

In this section, we work out the system density oper- 
ator under the action of the decoherence models of the 
previous section, and analyze how it leaves a character- 
istic trace in the cavity mode operators. 

A. Approximations 

We summarize first the approximations we make to 
arrive at an analytical solution. 
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(i) The most significant difference of our approach com- 
pared to related work on mobile mirrors is the "sudden 
approximation" : we assume that the duration r of the 
pulse is short compared to the system oscillation period. 
Since the damping rate re also determines the actual pulse 
length inside the cavity, we require 



Q <C min(l/r, re). 



(9) 



In this limit, the system's motion is "frozen" while the 
pulse is applied, and the Heisenberg equations for the 
photon mode operators can be solved without taking into 
account the dynamics of the system operators. 

(ii) At the same time, the pulse must be sufficiently 
long in order to restrict the cavity dynamics to the two 
degenerate modes mentioned previously. This is valid 
when the inverse pulse length is small compared to the 
cavity free spectral range v. Combining this with the 
"bad cavity limit" assumed below (rer 3> 1), but exclud- 
ing a too small cavity finesse, we have 



1/r < re < v. 



(10) 



(hi) Decoherence can be described by a dissipative 
functional as in Eq. JJ) if the system is weakly coupled 
to its environment and decoherence is happening slowly 
on the timescale set by the oscillation period 2n/Q,. This 
requires the inequality 

r th « n (ii) 

for the thermalization model (JHJ. For the pure spatial 
decoherence model (JJJ , we require that it takes more than 
one oscillation to increase the average system energy by 
one quantum hQ. This leads to 



where the rate 



r = 



D 



(12) 



(13) 



gives the depletion of the system ground state, see e.g. 
Folman et al. in [T^]. Given Eq. (|12(l . depletion happens 
slowly on the scale of the oscillation period. In both 
cases, the oscillator has a large quality factor. 



B. Decoherence 

We now show that the dissipative evolution can be in- 
tegrated in terms of the system covariances in position 
and momentum. For the equilibrium initial state we con- 
sider here, the mean values (x), (p) vanish at all times by 
symmetry. Using the spatial decoherence functional J7J, 
we get from the Liouville-von Neumann equation (JTJ 

d 



t:^ 2 ) = —(px + xp) 
at m 

— (p 2 ) = —mQ 2 (px + xp) 
at 



2D 



^-(xp + px) = — (p 2 ) - 2mVL 2 {x 2 ). 
at m 



(14a) 
(14b) 
(14c) 



Characteristic for spatial decoherence is that only the 
momentum width is increased by the diffusion coefficient 
D. This "squeezes" the system state in the phase plane, 
while the dynamics in the harmonic potential subse- 
quently leads to a rotation. The coupled equations (|14al - 
I14c|l can be solved, with the result (see, e.g., [l4j ) 



(x 2 ) t = (x 2 ) Q + 
(P 2 )t = {p 2 )o + 



hT 

2mVL' 
hmT 



- (2Qi-sin2Qi) (15a) 
(20* + sin 2Qt) (15b) 



(xp + px) t = — (1 — cos2f2i) , 



(15c) 



where the decoherence rate V has been defined in 11311 . 
We have used that equipartition holds in the initial state, 
{p 2 )o/m — mfl 2 (x 2 )o, which is obviously true for an ini- 
tial thermal state. Note that the decoherence Liouvil- 
lian does not describe a stationary solution in the limit 
t — > oo, this is because we neglected friction. 

For the thermalization model JHJ) , the variances can be 
computed similarly, and we find 



(x 2 ) t = (x 2 ) Q + 



E(Te) 



1 



-2r th t\ 



(16a) 



(P 2 )t = ip 2 ) +E{T e )m{l-e~ 2T ^) (16b) 
{xp + px) t = (16c) 



r 



in 



(16d) 



where E{T e ) = \Titt coth(S,0/2fc B T e ) is the average oscil- 
lator energy at equilibrium with the environment (tem- 
perature T e ), and where the rate T t h characterizes both 
the approach towards thermal equilibrium and the damp- 
ing of the system's average position and momentum. 

The key benefit of this formulation in terms of covari- 
ances is that it provides an exact solution for the system 
density operator |15L llq . The reasons for this are the 
thermal initial state we consider here and the Liouville 
functionals QIH) that are bilinear in x, p. We shall work 
with the Wigner representation W(x,p,t) of the density 
operator that has pro perties similar to a classical phase 
space distribution [lj|. For example, expectation values 
of symmetrized system operators S(x,p) are computed 
according to 



(S(x,p)) t 



dx dp 
2-Txh 



S(x,p)W(x,p;t). (17) 



Given the covariances, we find the Wigner function 

W(x,p;t) = Af(t)exp[-G(x,p;t)] (18a 

2G(x,p;t) 



7-^7 — h 7-57- (p — —7— ) (18b) 
{x 2 ) t (p 2 ) t V A t 



C t (xp + px) t 
A t ~ 2(x 2 ) t 



(x 2 )t(p 2 )t - \ (xp + px)t 



Af(t) 2 



(18c) 



(18d) 
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Note that for the normalization factor, one finds Af(t) < 
2 because of the uncertainty relations. 



C. Short probe pulse 

In a frame rotating at the cavity resonance frequency, 
the Heisenberg equations for the photon operators are 

d e = —ig (cos(2kx)a e + sin(2kx)a ) — na e + r/(t) + £ e 

(19a) 

a = ig (cos(2kx)a — sin(2fcir)a e ) — na + £ Q (19b) 

where £ e ,o are quantum noise operators that can be ne- 
glected as long as we calculate normally ordered quanti- 
ties, such as the intensity. In the sudden approximation, 
we assume that the system position operator x does not 
change during the pulse duration. It can then be treated 
as a constant that commutes with the photon operators. 

The initial condition just before time t p is vacuum for 
both cavity modes. The evolution of the odd mode is 
then given by (neglecting terms with vanishing expecta- 
tion value) 



a (t p + t) = —i sin(2kx) 

g — ge~ Kt cos gt — ne~ Kt sin gt 



g 2 + k 2 



V (20) 



for < t < t, where r is the pulse length, and the 
pulse shape was taken as a mesa-function: T)(t p + t) = r) 
for < t < t , otherwise it vanishes. With this choice, 
rj = \J k/(2t) din, where dj n is the boson operator of the 
probe pulse incident on the cavity. 

The average number of odd photons N Q at the detec- 
tor during the time window [t p ,t/] can be calculated by 
integrating 2na\{t)a {t) , see Eq. |JBJ. As mentioned in 
Sec. IIV Al we consider here the bad cavity limit where 
kt 3> 1. In this limit, we may set tf = t p + t, and the 
integration yields the simple result 



N o = R{a\ n a tn sin 2 (2fcr)> 



with 



R = 



g 2 K 2 



((J 2 



,2\2 



(21) 



(22) 



The average over the field and system operators in 
Eq. H21JI can be computed independently because we have 
assumed that the system+field density operator factor- 
izes at t — t p . The signal is thus proportional to the mean 
photon number of the probe pulse, (a\ n ai n ) = AT in . We 
discuss the signal fluctuations in Sec. IIV PI below. For 
the average over the system, the Wigner function l|18b(l 
gives after one elementary integration 



(sin 2 (2fcc)) tp 



dx dp 



sm 2 (2kx)W(x,p;t p ) 



2nh 

1 - exp [~8k 2 (x 2 ) tp ] 



(23) 



where the variance (x 2 )t has been calculated in Eqs. I|15a|) 
and Qlfiafl . 

We observe that the solution (|2(Jfl illustrates how our 
detection scheme conserves parity. Consider an incident 
pulse in a single photon state, \^>{t p )) — aj n |0), and as- 
sume that the odd detector clicks. This updates the sys- 
tem state to 



''odd click": <— > 1C sin(2kx)pA sin(2fca;) 



(24) 



where /C is the normalization. If the system has been in a 
state pa with definite parity, this state has the opposite 
one. Similarly, a click in the even mode detector updates 
the system state to a state with the same parity. In both 
cases, the "collapse" of the wave function does not lead to 
an a-symmetric state with less well defined parity. This 
ensures that symmetry of the system state can only be 
changed by decoherence. 

The factor R in Eq. ill'L'l) gives basically the probability 
that an incident photon actually interacts with the sys- 
tem. This can be seen from the number of even photons 
that is given by 



Af e = RN in + RN m (cos 2 (2kx)) 



(25) 



where R = K 4 /(g 2 + k 2 ) 2 . The first term is independent 
of the system position and gives the photons that did not 
interact with the system. This number reduces to 
for a vanishing coupling, g — > 0. The second term has 
a similar structure as Eq. 121|) . with the difference that 
here only the system operator with even parity occurs. It 
adds up with N Q to RNi ni similar to the two output ports 
of an interferometer. The fraction of 'useful' photons is 
thus R/{R + R) = g 2 /{g 2 + k 2 ). 

For typical experimental conditions that we discuss in 
Sec IV Al the probe wavelength is large compared to the 
width of the system position distribution ( "Lamb-Dicke 
limit"). We can then expand the exponential in Eq. i|23|) 
to get 



N ^4RN m k 2 (x 2 ) tp , 



(26) 



so that a larger signal is obtained with a shorter probe 
wavelength. This scaling breaks down, however, in the 
extreme case of A being comparable or shorter than the 
position width ("anti Lamb-Dicke limit"). This may also 
occur in the long-time limit, after heating has signifi- 
cantly broadened the position distribution. The signal 
then no longer increases linearly with (x 2 )t p - The expo- 
nential in Eq. vanishes, and the signal saturates at 
N Q = RNi n /2 regardless of the value of the (x 2 )t . In 
this limit, the probe pulse is no longer able to extract 
information about the system. 

We note that the same result N a = RNi n /2 may be 
arrived at on much shorter time scales (even for t p = 0) 
when the initially prepared system is already in the anti 
Lamb-Dicke limit. The system then heats as a result 
of photon scattering ("back action"). The exponential 



exp(- 



2 )o) in Eq. is in fact the Debye- Waller 
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factor for this scattering process. It gives the probabil- 
ity of the system still occupying the ground state after 
the photon pulse has impinged on it. In the anti Lamb- 
Dicke limit, the Debye- Waller factor is zero, and the sys- 
tem makes with probability unity a transition to a vi- 
brationally excited state. Since the detection of an odd 
photon is directly correlated, due to symmetry conserva- 
tion, to a transition from the even ground state to an odd 
state, it occurs for one half of all 'useful' photons RNi n 
that actually interacted with the system. 

One may also understand this result in the position ba- 
sis. In the anti Lamb-Dicke regime, the system position 
is smeared out over many probe wavelengths so that the 
phase of the backscattered light varies randomly from 
to 2ir. On average, one half of those photons that inter- 
acted with the system are sent into the odd detector. 



vibration frequency Q./2n 


50 kHz 


pulse length r 


20 ns 


cavity line width k/2-k 


50 MHz 


free spectral range v 


1GHz 


wavelength A 


795 nm 


Lamb-Dicke parameter (kao) 2 


0.073 


coupling strength g/2n 


10 kHz 


decoherence rate F/2n 


1kHz 



TABLE I: Realistic experimental parameters for a Rb atom 
fulfilling the validity conditions of the calculation, and yield- 
ing a signal to noise ratio larger than 3 in a spectral analysis 
(as discussed around Eq.(J2U). 



V. DISCUSSION 



D. Signal fluctuations 

For the discussion of the signal to noise ratio in 
Sec. IV Bl we also need the fluctuations of the odd de- 
tector signal l]2ip. 

The variance AN% can be computed from the factor- 
ized state pAF(t P ) as well. For a probe pulse in a coherent 
state, we find 

AN 2 = R 2 (N in (N in + l)(sin 4 (2^)) tp 

-N 2 n { S in 2 (2kx)) 2 tp ) (27) 

and 

{sm 4 (2kx)) tp = i (3-4e- 8fe2 <* 2 ^ +e- 32fc2 <*%) 

«3(4fcV) ip ) 2 . (28) 

In the last step, we made the long- wavelength expansion 
as in Eq. (|2fjj) . 

If the system were located at a fixed position, the terms 
proportional to N 2 n in Eq. (|27|l would cancel, leading to 
number fluctuations limited by shot noise. This is not 
the case here, however, because of the finite position un- 
certainty. For the Gaussian distribution at hand, the 
fluctuations of x 2 are of the same order and even some- 
what larger than their mean value (a: 2 ) itself. In the 
long- wavelength limit, Eq. I|27[) becomes 

AN 2 = R 2 (4k 2 (x 2 } tp ) 2 (2N 2 n + SN m ) ■ (29) 

For completeness, we also give the opposite limit of a 
short wavelength, although it is experimentally more 
challenging. Recall that the average signal N Q » RNi n /2 
arises because the phase 4> = 2kx is uniformly distributed 
between and 2tt, and the fraction of photons in the odd 
detector proportional to (sin 2 (j>) = | . The variance then 
becomes AiV 2 = R 2 {\N 2 n + §AT in ), where the second 
term is due to shot noise and the first is due to the vari- 
ance (sin 4 0} — (sin 2 (f)) 2 = A. This estimate agrees with 
the anti Lamb-Dicke limit of the general Eq. (|27|l . 



We show in this section that experimental conditions 
should exist where it is possible to spot the difference 
between spatial decoherence and thermalization by mea- 
suring the photon signals outside the cavity. 



A. Example I: Single atom 

The system we consider in the first example is a sin- 
gle rubidium atom trapped in a tightly confining mag- 
netic trap. The chosen parameters are listed in Table [I] 
The oscillation frequency is similar to those achieved 
with magnetic traps on atom chips |l2j |. For such trap- 
ping frequencies, ground state cooling leads to a spatial 
width (x 2 )l /2 of the order of a = {h/2mfl) 1 / 2 « 34 nm. 
Trapped ions cooled to the ground state with sideband 
cooling are also a good system for this purpose. We note 
that a cold thermal state would be sufficient as well. For 
the sake of this example, we choose a relatively long wave- 
length close to the Dl rubidium line. This gives a reso- 
nantly enhanced AC polarizability, while absorption and 
spontaneous emission can still be minimized with a de- 
tuning of several linewidths. The corresponding coupling 
strength g is estimated in the Appendix. Pulse length, 
cavity linewidth and free spectral range are chosen to 
comply with the approximations defined by Eqs. 10 and 

The last quantity, the decoherence rate T, depends on 
the specifics of the system. For an atom in a miniaturized 
electromagnetic trap on an atom chip, estimates for heat- 
ing due to magnetic field fluctuations are in the range of 
r ~ Is -1 12]. Other environmental perturbations can 
be added at will to enhance this rate in a controlled way. 
In fact, a tunable decoherence source is suitable for the 
unambiguous, experimental discrimination between de- 
coherence and thermalization-induced dynamics. 

The signal in the odd photon detector given by Eqs. 
(j2Tj) and is plotted in Figure M using the spatial 
decoherence model Q and the parameters given in Ta- 
ble The signal shows an overall increase because the 
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50 100 150 200 250 300 
tp [|is] 

FIG. 2: Average photon number N (t p ) hitting the odd de- 
tector (D) for Ni n — 10 9 incoming photons, as a function of 
the delay time t p in microsecond at which the measurement 
takes place. In the inset the initial oscillations are magnified. 
Parameters are given in Table (Rb atom) . The coupling 
to the environment is taken in the pure decoherence form of 
Eq.0. A thermalizing environment [Eq.JHJ] would not lead 
to the superimposed oscillations, but to an otherwise similar 
behaviour. 



environment heats the system and broadens its position 
distribution. The important feature of this plot are the 
superimposed oscillations. They stem from the breathing 
motion of the wave packet and are a telltale sign of spa- 
tial decoherence that affects position and momentum in 
a non-equivalent manner. Indeed, Eq. (|14bfl shows that 
decoherence only increases the momentum width which 
leads to a squeezed phase space distribution. The dy- 
namics in the harmonic potential makes this elongated 
distribution rotate at the frequency f2 so that its pro- 
jection onto the position or momentum axis oscillates in 
width at an angular frequency 2f2 (as expected from Eqs. 

Eli. 



B. Signal visibility 

Next, we discuss ways to extract the oscillations at 
2f2 of the odd detector signal that are characteristic for 
spatial decoherence. 

Let us consider probing times t p in the range where 
the system is in the Lamb-Dicke limit. That is, in Fig. El 
we are well before the saturated regime, i.e. t p < 100/is. 
The number of odd photons is then given by Eq. 12(j|) 
which, combined with Eq. (|15af> yields the peak-to-peak 
amplitude of the breathing oscillations 



8N m R(ka Q ) 2 -, 



(30) 



where (kao) 2 = Tik 2 /2mVi is the so-called Lamb-Dickc 
parameter. It is interesting to note that the oscilla- 
tion amplitude (|30[) depends neither on the initial system 
state nor on the delay time t p for the probe pulse [l8| . 



However, the last three factors, R, kao, and T/il are all 
less than unity, and a significant photon number can only 
be obtained with a bright probe pulse, N{ n 3> 1. In the 
example given above, Ni n ~ 10 9 for a 20 ns pulse (10~ 10 
J per pulse), and this large number compensates for the 
small fraction of 'useful photons' R » {g/^) 2 ~ 10~ 8 . 
The signal can be improved a lot with a larger coupling 
strength g. For the optimal value g = k, the ratio R 
takes its maximum value 1/4. The signal amplitude can 
also be increased with a shorter probe wavelength as long 
as one remains in the Lamb-Dicke limit. 

The contrast C osc of the breathing oscillations is de- 
fined by dividing the oscillation amplitude (|30l) to the 
smoothed background signal. Since that background in- 
creases with t p (in the Lamb-Dicke regime) , the contrast 
is time-dependent. In terms of the initial system temper- 
ature Tq, 



Cosc (tp ) 



2T/n 



coth{nn/2k B T ) + n r 



(31) 



where the dependences on photon number, detection 
probability and probe wavelength have canceled. The ini- 
tial contrast is maximized for low temperatures fcsTo < 
Ml when the system is cooled close to the ground state. 
Even then it is limited by the small ratio T/il, see 
Eq. i|12[) . For higher temperatures, the contrast decreases 
like HT /fceTo <C 1. The contrast reduction with increas- 
ing t p happens on the time scale coth(Tiil/2kBT )/T, 
where decoherence has approximately doubled the po- 
sition variance compared to its initial value. This time 
scale is of order 1 /T if the initial temperature is low. Oth- 
erwise, when the initial position distribution is already 
quite broad, it is of order fcsT /(?iiir), so the decoher- 
ence needs longer to double the initial width. 

To get a realistic estimate for the visibility of the oscil- 
lations, however, one has to take into account the fluctu- 
ations of the detector signal. The variance AN 2 is given 
in Eq. (|29|l for the Lamb-Dicke limit. It shows a super- 
Poissonian scaling with N 2 n to leading order so that we 
get a signal to noise ratio 



S/N(t p 



N 



Cosc (tp ) 



AA^o 



(2 + 3/7V m )V2 



(32) 



which is much smaller than unity. In order to resolve 
the oscillations, the S/N ratio has to be enhanced by 
repeating the experiment a number of times 



n cx « 10 /[S/N(t p )f 



(33) 



which ensures that the measured oscillation amplitude at 
an antinode exceeds the background noise level by three 
standard deviations. For short enough times this condi- 
tion yields a number of n cx ~ (fl/T) 2 3> 1 recordings for 
a given time t p . 

Alternatively, the oscillations can be extracted from 
the signal frequency spectrum in terms of the peak they 
give at 2f2. If we assume that the increase of the back- 
ground is negligible on the time scale 2tt/SI, the signal 
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fluctuations give an approximately white noise. In terms 
of the weight of the peak at 20, we find an improvement 
by a factor {T/ At) 1 / 2 with respect to the signal to noise 
ratio (|32|) . where T is the maximum probing time in the 
data set and At w t the time resolution. This shows that 
the sudden approximation (a small value of r) is actually 
required to extract the signal. We can take into account 
the slow increase of the background in an approximate 
way by writing the spectral signal to noise ratio as 

S/N (20.) « {T/T) 1 ' 2 S/N{t p = T/2) (34) 

(T/r) 1 / 2 2T/SI 

(2 + 3/N in y/ 2 coth(hQ/2k B T ) + TT/2 ' 

The signal to noise ratio shows a maximum for a specific 
observation time 



T, 



coth(hQ,/2k B T ) 



opt 



(35) 



which is T op t ~ 1/T for low initial temperatures. It is in- 
teresting that the system should be monitored up to the 
time where the contrast ll.il I) of the oscillations decreases. 
As a characteristics of the efficiency of the detection pro- 
cess, we obtain the maximum S/N 



max S/N(2ti) 



rtanh(MV2fc s T ) 



1/2 



(36) 



The spectral signal to noise ratio can be significant, i.e. 
we get max S/N(2Q) > 3 for the system parameters of 
Table 01 and it thus provides an evidence for the oscilla- 
tions. Note, however, that this signal is obtained from 
the records for the whole time series from to T opt . That 
is, a measurement repetition n cx w T opt /r w l/(rY) is 
necessary. On taking into account that a time resolution 
t < 0.1 r/fi 2 is required for getting such a significant 
S/N, we arrive at a necessary number of measurements 
n cyi « (£l/r) 2 3> 1. This condition is similar to the one 
we obtained for the resolution of the oscillations at a 
single point t p . This result also justifies the use of a con- 
trolled decoherence source to enhance the decoherence 
rate T. For example, with the parameters given in Table 
[fl about 10 4 measurements are required for evolutions in 
the range of to 100 fis to get the signal to noise ratio 
larger than 3. 



C. Example II: Nanoparticle 

Finally, we come back to the main motivation of this 
work and consider the decoherence of massive systems. 
As examples of such a nanoparticle we think of rece ntly 
realized nano electromechanical oscillators (NEMs) [H| 
or beads, which may be put in an harmonic confining 
potential by holding their edges or by suspending them 
in midair by electric (e.g. Paul trap in case they are 
charged), magnetic (in case they posses a magnetic mo- 
ment which maintains a meta stable state), or optical 
interaction (optical tweezers for small beads). 



The crucial parameter is the initial temperature which 
makes a difference with respect to the case of single 
atoms. While atoms can be routinely cooled down to 
the ground state of oscillation in a trap by various op- 
tical methods, there is no such efficient cooling scheme 
for arbitrary massive objects. We can safely assume that 
nanoparticles can be cooled down to below IK or even, 
in the near future, to the mK range, u sing e.g. special 
cryogenics or opto-mechanical feedback [2(j (for experi- 
ments, see Cohadon et al. in Ref. . 1 lj ) . With experimen- 
tally feasible trapping frequencies in the range of MHz, 
such temperatures are still large, Ml -C k B To, and the 
initial vibrational state of the system is highly excited. 

On the other hand, for large enough masses, such a 
temperature range is sufficiently low to reach the Lamb- 
Dicke regime, 



k 2 (x 2 ) 



2 2k B T 



«1, 



(37) 



where again (fccxo) 2 = Tik 2 /2mVl is the Lamb-Dicke pa- 
rameter (see Eq.ljSUJ). According to Eq. H23fl . this is 
the basic necessary condition for the proposed detection 
scheme to work. We can then use the general result of 
Eq. H36|) and expand the tanh function to lowest order: 



xaaxS/N(2Q) s 
and similarly from Eq. (|35(l : 



2k B T riT 



1/2 



opt 



2k B T 

nnr 



(38) 



(39) 



The requirements for the experiment can be quantita- 
tively estimated on the basis of these equations, regard- 
less of the specific implementation of the scheme. 

One important quantity is the number of measure- 
ments needed to get a statistically significant signal of the 
oscillations in the odd photon detector. First, a number 
of TLpt / r measurements has to be carried out to monitor 
the decoherence evolution from t p = to t p = T opt . Sec- 
ond, this full evolution has to be recorded a number of 
about 10/ (max S/N) 2 times to reach, by spectral anal- 
ysis, three standard deviation from the noise level. The 
total number of measurements is then 

n cx » 10/ (max S/ N) 2 x T opt /r ~ 10 {2k B T /hT) 2 , (40) 

which heavily depends on the initial temperature. It fol- 
lows that the decoherence rate should be artificially en- 
hanced to the maximum level allowed by the condition 

CU, i-e. r = o.i a 

Estimates for a reference nanoparticle are assembled 
in Table [H] For a trapping frequency SI = 27r 1 MHz, 
the decoherence rate can be as large as V — 2ir 100 
kHz. In order to maximize the signal to noise ratio, 
the pulse length can be as short as kt w 10. For a 
cavity linewidth k — 27r50 MHz, as the one in the Rb 
atom example, an initial temperature Tq — 1 mK yields 
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mass m 


< lcr 15 kg 


vibration frequency Q/2ir 


1MHz 


pulse length r 


30 ns 


decoherence rate T/2n 


100 kHz 


initial temperature To 


1 mK 


max S /N 


0.1 




100 fja 




10 6 



TABLE II: First block: parameters for a reference nanoparti- 
cle fulfilling the validity conditions of the calculation. Second 
block: calculated maximum of the 'signal to noise' ratio of the 
spectral analysis and the maximum recorded evolution time. 
Third block: the estimated number of experimental runs in 
order to achieve three standard deviations and the verification 
of the decoherence signal. 

the ratio max S / 'N (20) w 0.1, and we get a number 
?i cx ~ 10 6 where one measurement takes a time of about 
T ov t/2 ss 100/^s (not including the preparation time). 

In many possible physical realizations, the coupling 
to the nanoparticle is independent of its mass, e.g. the 
spring constant of the trap, K = mfl 2 , and the parameter 
D of the decoherence functional J7J). It is very interest- 
ing that T/Q = D / K and, hence, the maximum signal to 
noise ratio in Eq. I|38l) is independent of the mass of the 
particle. This proves that the scheme can be extended to 
monitor the decoherence of massive particles. The mass, 
in fact, can be scaled up to a limit m cr without degrad- 
ing the signal. This limit is related to our assumption 
that the particle remains smaller than the wavelength 
so that it couples to the light field via its polarizabil- 
ity. The upper mass limit is m cr = gX 3 ~ 10 _15 kg, i.e. 
about 10 10 Rb atoms, for a typical material density and 
visible light. The optimum time period to reach the max- 
imum spectral visibility, on the other hand, depends on 
the mass T opt ~ m2kBTo/ D and makes the necessary 
conditions of observing decoherence less demanding for 
smaller masses. 



VI. CONCLUSION 

To conclude, we have addressed the problem of observ- 
ing the decoherence process in massive objects. Exper- 
imentally quantifying the process for large masses and 
different environments is of paramount importance for 
the accurate theoretical modeling of this subtle transi- 
tion from quantum to classical. 

Specifically, we have shown that it should be possible 
to study the decoherence of systems that are trapped. In 
contrast to conventional spatial dephasing experiments, 
in which the observed object is freely propagating, the 
presented scheme will allow for more isolation, and hence 
for a better control over the coupling to the environ- 
ment. This is made possible by making use of a photon 
probe that scans the object for changes in its wave func- 



tion. Furthermore, the suggested experimental scheme, 
bypasses the need to create a clearly separated spatial 
superposition. This need presents an ever growing tech- 
nical difficulty for large masses in terms of the actual 
preparation procedure and in terms of the rate of deco- 
herence. A "pure" decoherence signal is apparent even 
when no superposition was created initially. The only re- 
quirement in the present scheme is to prepare the system 
in an equilibrium state of an harmonic potential. This 
is a further advantage over the ground state cooling re- 
quired by numerous other schemes. 

These features will greatly enhance the feasibility of a 
decoherence experiment with a scalable object mass. In 
addition, the interaction with optical cavity modes can 
be used to tailor the object's environment and to induce 
"decoherence on demand" [2l| ■ 

In this paper, we have limited ourselves to objects 
smaller than the wavelength that weakly scatter light. 
Future work will address moving mirrors that lead to a 
stronger optical signal. The fact that our scheme does 
not require strong reflectivity will allow the use of very 
light double sided mirrors (foils) with high transmittance. 
Preliminary results show that a similar Hamiltonian ac- 
counts for the coupling to the cavity modes so that most 
of the present analysis can be carried over, but with more 
favorable parameters. 

Finally, we have focused on environmental decoher- 
ence. However, if the system may be isolated well enough 
so that the coupling to the environment does not mask 
other weaker processes of localization, then one may 
perhaps be able to study also other proposed models 
|22L |23| . For estimates of the parameters required for 
opto-mechanical tests, see Bose et al. and Marshall et al. 
in Ref.0]- 
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APPENDIX A: INTERACTION HAMILTONIAN 

We describe the interaction between the electromag- 
netic field and the object by \aE 2 , where a is the po- 
larizability and E the (linearly polarized) electric field. 
In second quantization, keeping only the even and odd 
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cavity modes, the field can be written 



E 



huj c 



{a e cos(fcx) + a Q sm(kx) + h.c.}, (Al) 



where oj c is the (common) mode frequency, the wavenum- 
ber k — lu c /c, V is the cavity mode volume. In the rotat- 
ing wave approximation and adopting normal order, we 
hence find for the interaction 

xfko r 



H int = 77—77 { a t a e + a t a o + {a\a e - a\a ) cos(2kx)+ 

(A2) 



2e Q V 

+ (a\a + a\a e ) sin(2fca;)} 



We see here that the scattering from even into odd modes 
is accompanied by an anti-symmetric excitation of the 
object. If the photon leaves the cavity in the same mode, 
however, the object state is changed by the symmetric 
function cos(2kx). For the interaction Hamiltonian used 
in the main text, we have left out the part a\a e + aJa D 
involving the total photon number. Since this is a con- 
served quantity, it only contributes a global phase factor 
to the system+field wave function. 

Writing Tig for the prefactor in Eq. I|A2(1 . we get the 
coupling constant g. With an atom as system, we adopt 
a two-level model for the AC polarizability and cast the 
coupling in the form 



1A 



&TTUJA ^>A ~ 0J C A 



\ 2 



(A3) 



where lo a is the atomic resonance frequency andq I/7A 
the corresponding radiative lifetime. A is the cross sec- 



tion of the cavity modes and v = c/L the free spectral 
range. Eq. (|A3|I assumes an excitation not too far off 
resonance, \uja — u c \ <§; loa- Typical parameters are 
v ~ 10 9 /s, A/X\ ~ 10 3 , \oj a - w c |/7a ~ 10, and give 
g ~ 10 4 /s. This can be increased further with tighter 
focussing, smaller cavity lengths and working closer to 
resonance. A far off-resonant excitation with a shorter 
wavelength would improve the resolution of the position 
measurement. This is likely to be over-compensated by 
the smaller polarizability, since for u> well above a reso- 
nance transition, a oc 1/uj 2 — > 0. 

If the system is a nanoparticle with size smaller than 
the wavelength, the polarizability can be written in the 
Clausius-Mossotti form 



a((jj) = 47T£oO' 



£(w) + 2 : 



(A4) 



where where e(uj) is the particle's permittivity. The fac- 
tor fj, = (e(uj) — l)/(e(o;)+2) shows a resonance for exam- 
ple in metallic particles (collective plasma oscillation) . A 
value fx ~ 5 seems realistic while avoiding too large ab- 
sorption losses. This gives a coupling 



A* 
A 



(A5) 



where the 'critical mass' is defined in terms of the mass 
density g as m CI = pA 3 . Typical numbers (g — lg/cm 3 , 
A = 1 fim) give m cr ~ 10~ 15 kg. With the same numbers 
for the cavity mode as before, we get the fairly large value 
g k u(m/m CI ). 
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